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Wavefunction correlations and density matrices for few or many particles are derived from the 
properties of semiclassical energy Green functions. Universal features of fixed energy (microcanoni- 
cal) random wavefunction correlation functions appear which refiect the emergence of the canonical 
ensemble as A'' — > oo. This arises through a little known asymptotic limit of Bessel functions. 
Constraints due to symmetries, boundaries, and collisions between particles can be included. 
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. The standard tools of quantum chaos investigations include random matrix theory and periodic orbit theory 
(Gutzwiller trace formula), the Van Vleck-Morette-Gutzwillcr propagator, and many techniques and phenomena 
derived from these approaches. Standing somewhat to the side as an inspired insight is Berry's conjecture, which 

■ loosely stated is the idea that as — > eigenstates will be indistinguishable from superpositions of infinitely many 
(local) plane waves with random amplitude, direction, and phase, but with fixed wavelength appropriate to the local 

■ kinetic energy. In two dimensions, these assumptions result in strictly Gaussian statistics of the eigenfunctions and 
the autocorrelation function {ip* (x)ip{x + R)) = Jo(fca) where k is the local wavenumber and \R\ = a. 

The Berry random plane wave (RPW)[l| hypothesis is free of any specific dynamical information, except fixed 
^ ] total energy, which defines the "ensemble" (i.e. microcanonical). The perspective developed here suggests that by 
extending the RPW hypothesis we can conveniently accommodate many other constraints, incorporating information 
^H' about real systems. In fact this program has already begun, with Berry's inclusion of the presence of nearby hard 
walls [il, and Bies and Heller's soft boundary results [H, and multiple hard walls 3 . Related work by Urbina and 
RichterQ and one of us Q may also be viewed in this light. 

The idea of random waves subject to constraints is not confined to one particle in two dimensions. Indeed Berry 
gave the N - dimensional formula for free particles in his 1977 paperQ. Since the underlying idea in the RPW 
QQ [ hypothesis is uniform randomness within a quantum context, i.e. the underpinning of quantum statistical mechanics, 
00 • we must encounter some familiar territory as the RPW hypothesis is extended to the large N limit. In 1994, Srcdnicki 
^-H ' had suggested that the Berry random wave hypothesis was indeed a foundation for quantum statistical mechanics Q , 
■r^-j- , and showed that the appropriate canonical ensemble was reached for large N, depending on particle statistics. The 
' present paper shows more specifically what happens as the number of particles increases, through a nonstandard 
t , and apparently unpublished asymptotic form for Bcssel functions (we have not been able to find it in the literature, 
• although it "ought" to be there), which encodes the equivalence of the canonical and microcanonical ensembles of 
^ I statistical mechanics. In making the connections to quantum statistical mechanics one also needs procedures for 
incorporating constraints, which are an essential aspect of the theory. Thus our procedures for generalizing the RPW 
to include constraints, mentioned above, is an essential new feature, since the constrained eigenstates are no longer 
^ . random in Berry's (and Srednicki's) original sense. 
■ - ' Given a continuum at energy E, such as in an enclosure with walls very far away, we can perform the average over 
all random waves as a trace, i.e. 

inSMx')) ^Tt[5{E - H)\x){x'\], (1) 

which immediately yields Berry's result, apart from normalization which we choose differently here. However a trace 
over a basis is independent of any unitary transformation on that basis, so it does not matter whether we use a trace 
over a complete set of random waves or simple local plane waves; both give Jo(fca) for the case of one free particle in 
two dimensions. In this way the imaginary part of the retarded Green's function — ;^Im [G'^{E)] ~ 6{E — H) becomes 
central, formally convenient, and equivalent to Berry's RPW hypothesis. 
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II. PRELIMINARIES 



We begin by reviewing well known formalism to establish context and notation. The Green function completely 
characterizes a quantum system, whether it is interacting or not, or has few or many degrees of freedom. The retarded 
Green function G"*", i.e. 

G+ ^V^^-i^5{E-H), (2) 

where V stands for the principal value of the integral, is the basis for wavcfunction statistics and density matrix 
information, through the follow relations, with a convenient choice of normalization: 

< 7/.(x)^*(x') > = -hni{^\G+\^')/p{E) (3) 

= {^\6{E-H)\^)/p{E) (4) 

where 

p{E)=T^[5{E-H)] (5) 

and where < • • • > stands for the average over the degeneracies. We take these degeneracies to be of dimension up 
to ND — 1, where N is the number of particles and D the spatial dimension each particle lives in. (We use boldface 
notation, e.g. x for the N * D degrees of freedom.) If true degeneracies do not exist in a particular system, we 
can artificially open the system up to a continuum. For example, a two dimensional closed billiard does not have a 
degeneracy, but it acquires one if we open a hole in it and let it communicate with the outside unbounded 2D space. 
Of course this changes the billiard properties, and the size of the hole might be problematic, but in fact we shall 
never really have to open a system up in this way. The quantity 5{E — H) then implies the average over all scattering 
wavefunctions at fixed energy E. 

There are other interpretations which can be put on the average correlation < '!/;(x)?/'*(x') >; for example we can 
imagine a large number of potentals which differ in some far away place, and in a way so as to all have an eigenvalue at 
a particular energy. Then, the average has the interpretation of the average over this "disorder" ensemble. A slightly 
different procedure is advocated by Richter et. at, wherein an energy average is takenQ. Another interpretation can 
be applied to individual eigenstates in a closed system, assuming they are at least locally uniform in their properties, 
by taking the average over different points of origin x. This is particularly appropriate when the analogous classical 
system is chaotic, as mentioned above We will be evaluating the Green functions scmiclassically in what follows, 
restricting the time over which the contributing trajectories propagate. 

The wavcfunction correlation is equal to the coordinate space matrix element of the constant energy density matrix: 

< ^(x)^*(x') (x|J(i;-i/)|x')M^) -p(x,x',S) (6) 
Reduced density matrices can also be derived from wavcfunction correlations ; e.g. 

p{xi,x[,E)^ dx2dx3---dxN p{xi,X2,--- ]x[,X2,--- -jE), (7) 



the one particle reduced density matrix. 

We can approach the correlations via Fourier transform from the time domain, since 

oo 

S(E-H) = — f e'^'/'^e-'"'/'' dt. (8) 
27rAi J 

— oo 

Thus the statistics, density matrices and correlations are derivable without further averaging by knowing the time 
propagator. 

In the following, we define the Green function propagator G'(x, x',t) and the retarded Green function propagator 
G+(x,x',i) as 

G(x,x',f) = (xle-^^^/'^lx') 
G+(x,x',i) = ^eW(x|e-^^*/''|x') (9) 
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where 0(<) is the Heavyside step function Q(t) = 0, t < 0, Q{t) = 1, t > 0. It is very rewarding to expand the 
propagator in semiclassical terms, involving short time (zero length) and longer trajectories. We take Gdirecti^,^ + 
r, t) = (x| exp[—iHt/h]\x+r), the very short time semiclassical propagator, which for N particles each in D dimensions 
reads 



G 



(x, X + r, t) « f e™-V2nt-»v(x+i)t/;. (^q) 



direct 



where = |rp. 

It is not difficult to cast the Fourier transform of this short time version to fit the definition of a Hankcl function, 
i.e. 



where d ND/2—1, k = fc(x+r/2, E) and H^]^\kr) — Jd{kr)+iNd{kr) is the Hankel function of order d, and Jd is the 
regular Bessel function of order d. The wavevector k varies with the local potential, i.e. h'^k{x, EY /2m = E — V{'x). 
Here, using only the extreme short time version of the propagator, wc must suppose r is not large compared to 
significant changes in the potential, but this restriction can be removed by using the full semiclassical propagator 
rather than the short time version. For the case of one particle in two dimensions, d = 0, and we recover Berry's 
original result for one particle in 2D, {tp* {x)tp{x + f)) oc Jo(fcr). 
According to the short time approximation, for any N, 

where k — k{x,E). This result includes intcrparticlc correlations through the potential V{x.) and the spatial depen- 
dence of fc = fc(x, E); the diagonal r = limit (following section) is equivalent to classical statistical mechanics. The 
implications of this for the nondiagonal short time Green's function are intriguing. The way r is defined, it does not 
matter whether one particle is off diagonal (x^ ^ x^/) or several or all of them. For given r, the Green's function will 
be the same, apart from changes in the potential y(x + r/2). 

It is interesting that although the short time Green function is manifestly semiclassical, the energy form, e.g. Eq.fT^ 
is obtained by exact Fourier transform of the semiclassical propagator, rather than by stationary phase. 

III. DIAGONAL LIMIT 



The diagonal (r — > 0) N body Green function is obtained using the asymptotic form 

\im. Jd{kr) = ^ 
r^O I (a + 1) 

we obtain 




U"! -hikr) = I I Ti^ I ^ I (13) 



where the second form uses Stirling's approximation, n\ ~ n"e~"-\/27rn, and is appropriate below when we consider 
large N . Wc note that this behaves as k^'^ ^ {E—V{x)Y. This factor is familiar from the computation of the classical 
density of states. Tracing over all x results in 

ra 1 f ^^^^ f dx.dp 



27rft2 T{d+1) yAn 



J '^5{E-Ha{p,x))=pa{E) (15) 



i.e. the classical density of states. The association of the short time propagator with the classical Hamiltonian and 
classical density of states is well known. The Berry RPW hypothesis, the short time propagator, and the classical or 
Wcyl (sometimes called Thomas-Fermi) term in the quantum density of states are all closely related. 

The quantum spacial integral is over all coordinates, so how does the classical partition function emerge if the 
classical integral is only over classically allowed coordinates? For forbidden positions, k is imaginary and can be 
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written as say in. An identity for Hankel functions can then be used {i'^^^ Hn \ix) = ^Kn{x)) to show that the 
green function is real so that the imaginary part is zero, explaining why the integral is only over classically allowed 
positions. 

As long as r = (i.e. diagonal Green's function) the results obtained within the short time propagator approximation 
for any quantity in the presence of a potential {including interparticle potentials such as atom-atom interactions) will 
be purely classical. Since we will be discussing the equivalence of the results from the different ensembles for r ^ 0, 
it is useful to recall how the classical coordinate space densities in the different ensembles can be shown to coincide 
since this corresponds to the r = case. 

The normalized phase space density in the microcanonical ensemble and the phase space density in the canonical 
ensemble are given by 

p,i{p,^,E) = ^—S{E-Ha{p,^)) (16) 
Pci{E) 

and 

p,,(p,x,/3) = -i—e-^^-(P^-) (17) 

Qci(p) 

respectively. The density of states and partition function are of course the normalization factors so that 

Pci{E) = J dicdp 5{E ^ H,i{p,x.)) (18) 

QciiP) = J rfxdp e^^^-^P-'') (19) 
Integrating each phase space density over momentum space allows us to compare the coordinate space densities: 

Pe.(x,i^) = (20) 
^-'(^'/^) = Jrf^e-mx) (21) 



with p = ^2m{E- ^(x)). 

Using the relationship between E and (3, E — {V) ~ ^j^, where (V) is the ensemble average of the potential in one 
of the statistical ensembles, the coordinate space density becomes 

^ (2..(d+l)//3)-^(l+ (<^^^-|j-»^ )' (22) 



In the limit N oo { d ^ oo) this is 



= (2m(d+ l)//3)''e«^>-^(^»'' (23) 

(24) 



p2d p~V{x)i3 



J dx p2d J e-^(^')P 

This is one of the standard ways of establishing a connection between the ensembles Q. 

Since the diagonal Green's function gives classical results we can use it to study classical properties. For example, 
we can inquire about the average two particle spacing distribution pE{f 12) or the probability density for a single 
particle Pe {xi ) starting with the short time semiclassical Green's function and the results will coincide with classical 
microcanonical statistical mechanics. This statement holds for all N . Similarly, in the large N limit the canonical 
ensemble results for these quantities must emerge. This point becomes more interesting for the non-diagonal case, 
considered next. 

IV. LINK TO THE CANONICAL ENSEMBLE 
A. Bessel functions become Gaussians 



As yet we have found nothing too surprising or useful beyond standard classical statistical mechanics. This changes 
when we consider the large N limit for the non-diagonal Green's function, r 7^ 0. Taking the large N limit of Eq. [T^l 
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FIG. 1: As N increases, the combination -^Jdix), where d = ND/2 — 1, approaches a Gaussian. This is the key link between 
the quantum microcanonical and canonical ensembles. 



we are confronted with a new question about Bcssel functions. The large d limit of Jd{x) is indeed well known, but 
this is not yet sufhcient for our purposes. It reads 

V -JIZXT = o,.T^^, , 1^ ~ ^ (25) 



This is the standard formula given in the usual references. Eg. 1251 should be the first term in a power seres for Jd{kr) 
in kr. Another standard result is the power series expansion, valid for all d and kr: 

°° (-1)™ / kr^ 2m+d 



m=0 



We actually require a different asymptotic result. What make our demands unusual is that, assuming we want the 
energy to increase in proportion to the number of particles (appropriate to many applications of the large N limit), 
then k ^ Ve ~ ^/N ^ Vd; this means that for fixed r the combination (kr) is increasing as Vd as d — > oo. If the 
argument of the Besscl function increases without bound along with it's order, some new considerations come into 
play. We find the desired form using Eg. 1261 after summing a series recognized as that of a Gaussian Taylor expansion, 

where again /i^fc^/2m = E — V^(x). Note that as d ^ oo, the argument of the Gaussian holds fixed because of the 
factor of d + 1 in the denominator of that argument. Figure [T] illustrates the convergence to the Gaussian as N 
increases. The asymptotic limit in Equation [27] is not in the usual references, although related results have been given 
for N-bead polymer random chain end-to-end distributions [loj. The connection between the path integral for the 
propagator and polymer chains is well known 

It is interesting that a Gaussian emerges from Bessel functions in the large TV limit. We can put Eg. \T7\ together 
with Eg. [T2land Eq. [H and express the result, as ^ oo. 

For noninteracting particles moving in zero potential but confined to volume V the short time approximation becomes 
exact and k is constant. For this system the wavefunction correlation becomes 

< V^(x)V'*(x + r) > = p(x,x',S) ^ J_ ^-k\VMd+i) _ (29) 
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Something familiar is emerging, here derived in the unfamiUar context of fixed energy (microcanonical ensemble) . For 
comparison we recall the standard result for the ideal gas at temperature t[T^: 

Mf^;^. = J^e-.V.= ,30, 



where A = hj ^phrmnT is the thermal wavelength. Indeed for the free particle case, k is fixed by E and (/v) = 
Dl2NnT = h^k'^/2m, where K is the kinetic energy and k is Boltzmann's constant. 



The canonical ensemble result for the propagator has "dropped out" of the asymptotic large N limit of a micro- 
canonical Green function, at least for noninteracting particles, and an unusual asymptotic form for the Bessel function 
has emerged as the link. With some caveats, the statement 

S{E-H)^e-(^" (32) 

has meaning in the large N limit, where it is understood E grows as N, and a temperature extracted. At a qualitative 
level, Eq. [35] merely expresses the known equivalence of the ensembles. In the case of an interaction potential, the 
relation between E and temperature is of course problematical. 



B. Interacting Particles - Short Time Limit 

We can say more about interacting particles using only the short time propagator introduced above. Longer time 
events will be discussed in Sec. I VII The short-time approximation to the correlation function for large N, which is 
equal to the matrix elements of the density operator in coordinate space using our normalization, (Eq. I28p is given by 

PcKx,x',i^) . ( e-^^V4(..i) (33) 



with hk = 2m{E — V{^^^y^)) and r = |x — x'|. Again, the Gaussian form of this expression arises from the 
asymptotic limit of the Bessel function. In the interacting case this can again be brought into the same form as the 
equivalent expression at constant temperature: 

In order to make the connection we must identify the energy with a certain temperature. This relationship between 
E and (} is 

ND 

B-iV)^^ (35) 

where {V) is the ensemble average of the potential in one of the statistical ensembles. Using this relationship in Eq.[33l 
gives 

1 m f k^V _^uv)- 



p.(x,x',i^) = -^^^-j e-^e--^ (36) 

In order for Eq.[3n]to be equivalent to Eq.[M]the term with {V) — V must be negligible. This is true for configurations 
of particles which possess the typical (and vastly most probable) sum total kinetic energy for all the particles. Since 
the typical total kinetic energy is by far the most probable, nearly all points in configuration space lead to small 
values of {V) — V , and that term is negligible almost always. The remaining terms in Eq. [36] and Eq. [34] are shown 
to be the same by the equivalence of the classical ensembles as shown in Sec. IIIII 

It is also telling to trace over the coordinates of all but one of the interacting particles, given by a coordinate y. 
We thus seek the reduced density matrix, diagonal or off diagonal in y. The trace will over many coordinates be 
overwhelmingly dominated (in the large N limit) by the most probable total kinetic energy for all the particles. Then 
we find 

G{y, if, (i) - X-3N-2^-n^-/x' (37) 



where — \y — t/^ and A = h / \/ 2'KmKT . Thus the quantum mechanical single particle Green function and density 
matrix make sense as their imaginary time counterparts in the N oo limit, in accordance with well known results 
for the canonical ensemble. 
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C. Large N limit and Boltzmann averaged Green functions 

Even though it is a necessary consequence of the equivalence of the ensembles, it is interesting to establish the 
generality of the Boltzmann average over the energy of a noninteracting subsystem in the following way. Suppose 
N — M particles are no longer interacting with the remaining M particles, but their states are correlated by having 
been in contact in the past with the total energy fixed at E. In the time domain and in an obvious notation we have 

G+ (y, z; y', z', t) = ih G+ (y, y', t)G+ (z, z', t) (38) 

Then the Fourier convolution theorem can be apphed to the Fourier transform into the energy domain, i.e. 

oo 

G+(y,z;y',z',i?) = ^ j G+_,,iy,y',E~E')G+jiz,z',E')dE' (39) 

— oo 

which incidentally leads to some rather unlikely looking identities for Bessel functions; the reader may easily generate 
them. Our purpose is served if, focussing on the subsystem of M particles, we trace over the N — M y coordinates. 
This gives 

times a volume factor, in the case of an ideal gas. The second term is not a function of E'. Therefore the integral of 
it times Ga/(z, z', E) is proportional to d{z' — z). So long as z 7^ z' that term is zero. Neglecting all unimportant (for 
this argument) factors this leaves 

Tiy [G _ J^J {E - E' )] {E - E' f ^E'^'^-''' i I- — ] ^ g-/3is' (4;^) 

with of course /3 = 1/kT. In arriving at Eq. 1411 we used E = §NkT for the case of particles embedded in D 
dimensions. Finally we arrive at 

00 

Try[G+(i?)] (X J e-^""' G+,{z,z\E') dE' ^G+,{z,z',l3) (42) 
—00 

in the large N limit. This establishes the generality of the Boltzmann average over the subsystem energy for large N. 
This discussion establishes again the connection between the canonical and microcanonical ensembles, however in a 
way not involving the Bessel functions and their asymptotic form, so it is less general than other results in this paper 
valid for any N. 



D. Stationary phase canonical limit 

It is also possible to recover the Gaussian form in Eq.[2H]by carrying out the integral in Eq. [TT]by stationary phase, 
provided the real factor involving t in the denominator is taken into the exponent, as —ND/2\ogt i.e. 

00 

G+(x,X + r,S) = — / ^J!l-y^\"nr-mt-^yi.+ ^)t/n+^Et/n-ND/2logt (43) 



The complex stationary phase point t* in the large N limit becomes t* = —iNDh/{2{E — V)) , yielding the same result 
as in Eq. 1281 with h?k(x., E^ /2m = E — V{x.), and making this another route between the quantum microcanonical 
and canonical ensembles. Since the positions are arbitrary we cannot however identify the average kinetic energy with 
E — V, and thus without further averaging we cannot associate t* with any inverse temperature. It is interesting 
nonetheless that there is a complex time t* appropriate to every position x, even if that time is not related to the 
temperature. For an ideal gas the stationary phase time is t* ~ —ih/nT = —i(3h, after making the identification 
E = ND/2kT. A discussion about traces over most of the coordinates and the recovery of the usual temperature 
through {K) = D/2NkT proceeds as in Sec. HVB] 
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V. CONSTRAINTS 

111 the large N limit the ergodic hypothesis is strongly motivated, but statistical mechanics does not pre-suppose 
that ergodicity is unchecked; rather constraints are always present, such as walls and boundaries which control volume. 
Ergodicity is then defined with respect to these constraints. The guiding idea in this paper, i.e. the extended Berry 
RPW hypothesis, is that eigenstates of the full system are "as random as possible, subject to prior constraints" . In 
this way thermodynamic constraints arise naturally. The real time, real energy (microcanonical ) semiclassical Green 
function approach not only automatically generates the averages required to get appropriate wavefunction statistics, 
it also provides a natural way to include many constraints such as walls, symmetries, and even the existence of 
collisions between particles by going beyond the short time limit term to include returning (not necessarily periodic) 
trajectories. The semiclassical Ansatz for these extended problems in the presence of constraints is 

G(X, X', i) « Gd^reM, x', i) + ^ G, (x, x', t) (44) 

j 

where Gj(x, x + r,t) is a semiclassical (Van Vleck-Morette-Gutzwiller) Green function, 

expfz5,(x,x';t)/n-^j (45) 

corresponding to the j*'' trajectory contributing to the path from x to x + r, and Gdirect{^, x + r, t) is given by Eg. 1101 
The short time term G direct (x, x + r, i) , is singled out as the shortest contributing trajectory: supposing r to be small 
compared to distances to walls etc., we still have a short time, ballistic trajectory as quite distinct from trajectories 
which have traveled some distance away and come back. There are cases where this separation is not clean; for such 
cases we can adjust notation accordingly. Note that since a trace over all position is not being taken, there is no 
appearance semiclassically of periodic orbits as the only surviving contributors. "Closed" orbits however can play a 
large role semiclassically, a fact recognized long ago by Delosfisj. 

A. A'^ particles and a wall 

A very useful example is provided by a plane Dirichlct wall felt by all the particles (e.g. i/'(a?i,X2, ■ ■ -xn) — for 
Tji ~ 0,1 ~ 1,- • -N), as in a gas confined by a rigid container. The Green function and eigenfunctions must vanish 
if one or more particles approaches this wall. We can use the method of images, generalized to N particles, if the 
particles are noninteracting. (The interacting case can in principle be handled by semiclassical trajectory techniques 
which we bring up in the next section.) 

The Green function Guja/i(x,x') will consist of the shortest distance contribution for which all particles take a direct 
path from x to x', plus paths where one particle has bounced off the wall, paths where two particles have, etc. These 
histories are included automatically if we apply the symmetrization operator which imposes the image reflections. 
This operator can be written 

N 

i i i<j 

where Ri is the operator for reflection about the y = axis for the i*'' particle. Applied to the Green function 
G(x, X + r, i), considered as a function of the coordinates in x in the absence of the wall, TZ yields the series 

Gn^all (X, X' , t) = G direct (x, x' , t ) - ^ G, (x, x' , t) + ^ G„- (x, x' , i ) (47) 

i i<j 

where Gi(x, x',t) corresponds to the i*^ particle getting from Xi to x'^ by bouncing off the wall while the others 
take direct paths, etc. The Fourier transform gives an analogous equation for Gwaiii^,^' , E). The effect of the 
symmetrization is to create Green function sources reflected across the wall and given proper sign, in the manner 
familiar from the method of images. The short time path is shown by the direct path solid line in Fig[2l corresponding 
to the term Gst(x, x',t). The bounce path is equivalent to a source reflected across the wall with an opposite sign, 
i.e. the method of images. Define 

-il„[Gl(x,x + r,i5)] .^(i!)" ^^^mFJM) (48) 



Gj(x,x';() = 



27Tih 



I\ ±J / ^ 



Det 



9x9x' 



9 




FIG. 2: A short and a bouncing path for a particle propagating near a wall. The bounce contribution, if viewed by the image 
method, is equivalent to a contribution of opposite sign coming from the reflected point with the wall removed. 



Then 



(49) 

This is the general result for any N . It would appear to be difficult to take it further, since all the distances, e.g. 



where is the reflected particle coordinates, involve square roots. However if we use the large N asymptotic 
form, we find, using Fd{kr) exp[— fc^r^/4(ci + l)]/2'^d!, 



where 7 — k'^ /A{d + I) — tt/A^ and A| — [rfY ~ "^1- Since is the "direct" distance from Xi to x-, (see Fig [2]), 
records the distance change upon reflection of the particle. We note that (and thus the Green function) 
vanishes as any particle approaches a wall in either x or x'. It is also simple to see that the single particle density 
p{x) in this noninteracting case becomes, for large iV, 

p(f)=Po(l-e-4^-') (52) 

where x is the distance to the wall and po is the density far from the wall. 

The formulas Eq. 051 and Eq. [HI] generalize Berry's result for the wavefunction squared of one particle in two 
dimensions near a wall, namely 

The Gaussian we get for large N has a very simple interpretation. First we note that for noninteracting systems in 
the canonical ensemble we can write the total density matrix as a product of one particle density matrices. This is 
essentially the form of Eq. 1511 since we can write each one particle density matrix as 

f 1 _ ^~l{\3''-x'\^-\3-3'\^)\ ( 1 _ ^-^\3^-S\^\ 

p{x!,x',P) = e-^l-^-^^'l^/^^^— ^ - -4-- —WzU (54) 



Jdx{l- e-7l5«-sP) Jdx{l- e-7|5"-sP) 
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FIG. 3: The particle symmetry or antisymmetry condition is equivalent to requiring mirror symmetry or antisymmetry across 
the Xi — Xj (hyper)plane. This corresponds to having additional contributions from the images of the particles reflected over 
the symmetry planes. 

where the second form is the diagonal element. However Eg. 1541 also arises as the density matrix obtained from the 
Boltzmann average of Berry's result; i.e. averaging the fixed energy results over a canonical distribution of energies, 
as can be seen from the integral 



Jk {1 - Jo{k\x^ ~ x\)) e-'3'*''='/2'"dfc 

_0 

oo 





1 - 



(55) 



For D = 2 and = 1 a Boltzmann average yields the Gaussian. Indeed this necessarily holds in any number 
of dimensions; i.e. the appropriate Boltzmann average of J d{kr) / [kr)'^ must yield a Gaussian for any d. In the 
thermodynamic N oo limit for noninteracting particles, each particle separately is Boltzmann distributed over 
energy, so the result must be the same as a Boltzmann average of the one particle results for any dimension D and 
for any constraints. 



B. Symmetries - Fermions and Bosons 



Particle symmetry is an essential part of the many body problem. It's effect, like other symmetries, is to generate 
permutations where the distances have changed due to particle exchange. Figure [3] shows this effect graphically. It is 
gratifying to see directly that permutations which induce large new distances (coming from remote pairs of particles, 
where "remote" is a relative term depending on the temperature) make little contribution. Consider N noninteracting 
Fermions or Bosoiis; we wish to compute the reduced density matrix for two Fermions or Bosons. This is a well known 
result for TV — > oo|12l|. The symmetric or antisymmetric Green function is 

where r„ = \J\xi — Xp/p + • • • + |a;]v — ^pn'Vi {Pii ' ' ' tPn} is the nth permutation of {1, • • • , N}, and e„ = 1 if the 
parity of the permutation is even and e„ = ±1 if the parity of the permutation is odd (with the upper sign for bosons 
and the lower sign for fermions). 
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(V' [xi - ■ -XNjVKXi - ■ -xn)) = TT^T 



TV! / l2 \ 



In the limit that N is large, this becomes 

1 TO /fc2\ 

The diagonal component of this with the r„'s written out explicitly is 

/ ,2\d N\ 

(r(.l...xlVWxl....lV)) = (^) E^" (59) 

Up to the normalization constant this is the constant temperature density matrix for N noninteracting fermions or 
bosons: 

Again the identification E = nT was used. This can be rewritten as an integral over wavevectors: 

(|V'(x)|2) = AEe« / dki---dkN e-^'"'''=i'/2™+^'=i-(^i-^pi) • . . e-'^''''="'/2m+*fc„-(5iv-5piv) (gi) 



where A = 2 (E)N\7:h:^d\ \Tki ( ^rm ) normalization constant. Rearranging gives 



TV! 



If the volume that the particles are confined to is large but finite, 

JV! 



{\ll;{x)\^)dx3...dxN = AV'^-^ J^Cn dii e-P^ ^ /2rn^^{k^-k^^)■S^^^{k2~k^2)■S2g^^ (gg) 

n •' 

For fermions if the wavevector of any two particles are the same the term is killed by the term with the wavevectors 
reversed in accordance with the Pauli principle. This leaves only two terms 

r N\ .... 

{\iP{^)\'^)dx3---dxN = Ay^-^E^" / dk e-''''''''/2"e*('=i-'=''i)-^ie*('=^-'=''^)-^^ (64) 

For bosons there are also only two types of terms, but each is multiplied by the same factor since like terms are added 
together. Either way, carrying out the integral over k, 

J ^I'^WI^^'^^^^ • • • ''^ - Jdx:dx.il±e--i-^-Jfn ^''^ 
This is the well known result for the density of two noninteracting fermions or bosons. 
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VI. SCATTERING 



A hard wall is a potential energy feature which induces a boundary condition, requiring the wavefunction or Green 
function to vanish as the wall is approached. Softer potentials do not induce fixed boundary conditions and require 
a different treatment. A potential may still however be thought of as a constraint: we consider waves as random as 
possible subject to the existence of a potential, be it fixed or interparticle. In practice this means we return to the 
Green function formulation used throughout. 

Consider a soft repulsive or attractive potential somewhere in a noninteracting gas. Assuming no boundaries, 
mutually noninteracting particles can interact with the potential or 1 times. (We assume for simplicity that the 
potential is short ranged. Because of the ergodicity assumption inherent to the random wave hypothesis, the presence 
of remote walls would actually make no difference.) This circumstance develops along lines very similar to the wall, 
except that we cannot use the method of images. It illustrates the use of the full semiclassical propagator within this 
formalism. 

Eq. |46] and Eq. [47] both hold, with the effect of Ri changed to mean "the particle takes the path from initial 
to final coordinates in which it deflects from the potential, if such a path exists classically" . For N particles, there is 
a "direct" term in Eq. 27] where no particle interacts with the potential, N terms where one of them does, etc. We 
have, in the simple case shown in Fig. 2] and in analogy with Eq. ITfl 

G(x, x', i) — Gdirecti^j ^' jt) + y ] Giiounce.ii^j ^'j ^) + ^ ] Gbounce,i,j (x, x', t) + • • • (66) 

with Gdirecti^,^' ,t) given by Eq. (TUl and e.g. 

(W-l)D / -, \ a2Q f ' -f-^ 5 

Cx6o™(z,y„z + r,y„tj~ ^ J y^^^^J e [bt) 

Considering this term where only the i*^ particle with coordinate interacts with the potential, we have — 1 
"spectator" z particles, and the propagator becomes a product of the noninteracting Green function for A^ — 1 
particles and a more complicated Van Vleck semiclassical term for the colliding particle. The noninteracting part 
contributes a term (A^ — l)D/2\ogt in the exponent along with the one particle classical action of the i*'' particle. 
For sufficiently large A^ , and tracing over the z particles, this factor leads again to the usual time condition t* = —i/3h 
and a thermal average of the one particle energy Green function under the Fourier transform from time to energy, as 
in Equation 221 

G(y, y', S) « G(y, y', /3) = Gd^rectiy, y', /?) + ^ Gbo«„ce,,:(y, y', /3) + ^ Gbo«„ce,»j (y, y', /?) + •• • (68) 
t* — —i(3h becomes the imaginary time over which the action for the y coordinates are evaluated. 



VII. CONCLUSION 



Starting with Berry's random plane wave conjecture for chaotic Hamiltonian systems, we have followed it's impli- 
cations for moderate and large numbers of particles A^. In the large A^ limit we have necessarily arrived at some 
familiar territory in statistical mechanics. We have adopted a Green function, semiclassical perspective, arriving at a 
Gaussian-Bessel function asymptotic result for energy Green functions, providing an analytic connection between the 
quantum microcanonical and canonical ensembles. We have extended the incorporation of constraints into the random 
wave hypothesis, considering several types of constraints, including walls and interparticle collisions. Indeed the guid- 
ing perspective has been to make quantum waves "as random as possible subject to known prior constraints" . This 
must ultimately be equivalent to the ergodic hypothesis of quantum statistical mechanics. The nonstandard methods 
and perspective used here may possibly lead to new avenues of inquiry, and it is our hope that the semiclassical 
approach might permit new ways of treating strongly interacting systems. 

The next stage in the development of this approach is to consider short ranged potentials between particles, i.e. 
interparticle collisions. The first corrections to the free particle limit involve binary collisions, which can be computed 
semiclassically or using a delta potential appropriate to s-wave scatterers. Again the effect of the other particles will 
be to provide a thermal reservoir which essentially averages the Green function over a thermal distribution of energies 
(if N is sufficiently large). We save this for a future paper, where we hope to examine specific potentials and derive 
two particle radial distribution functions. 
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FIG. 4: A short ballistic and a colliding path both lead to the same final point for a particle propagating near a localized 
repulsive potential. The colliding path cannot be treated by the short time approximation; rather, a Van Vleck Green function 
is required. In this term, all but the i*** particle remain in place. 



Acknowledgments We thank Adam Wasserman for helpful discussions, and the National Science Foundation under 
grant NSF-CHE- 0073544. 



[1] Berry M V 1983 in Chaotic Behaviour of Deterministic Systems ed G looss G, R Helleman and R Stora (New York: 

North-Holland) p 171 
[2] Berry M V 2002 J. Phys. A: Math. Gen. 35 3025 
[3] Bies W E and Heller E J 2002 J. Phys. A: Math. Gen. 35 5673 
[4] Bies W E, Lepore N, and Heller E J 2003 J. Phys. A: Math. Gen. 36 1605 
[5] Urbina J D and Richter K 2003 J. Phys. A: Math. Gen. 36 L495 
[6] Heller E J 2006 Mol. Phys. 104 1207 
[7] Berry M V 1977 J. Phys. A: Math Gen. 10 2083 
[8] Srednicki M 1994 Phys. Rev. E 50 888 

[9] Jancel R 1969 Foundations of Classical and Quantum Statistical Mechanics (Oxford: Pergamon) 
[10] Kleinert H M 1995 Path Integrals in Quantum Mechanics, Statistics, and Polymer Physics 2nd edition (Singapore: World 
Scientific) 

[11] Chandler D and Wolynes P G 1981J. Chem. Phys. 74 4078 

[12] Pathria R K 1996 Statistical Mechanics 2nd edition (Oxford: Butterworth-Heineman) 
[13] Wang D M and Delos J B 2001 Phys. Rev. A 63 043409 
[14] Berry M V 2002 J. Phys. A 35 3025 



